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Abstract. We describe an approach to onhne analysis of micro-rheology data using a 
multi-scale time-correlation method. The method is particularly suited to process high- 
volume data streams and compress the relevant information in real time. Using this, 
we can obtain complex moduli of visco-elastic media without suffering from the high- 
frequency artefacts that are associated with the truncation errors in the most widely 
used versions of micro-rheology. Moreover, the present approach obviates the need to 
choose the time interval for data acquisition beforehand. We test our approach first 
on an artificial data set and then on experimental data obtained both for an optically 
trapped colloidal probe in water and a similar probe in poly-ethylene glycol solutions at 
various concentrations. In all cases, we obtain good agreement with the bulk rheology 
data in the region of overlap. We compare our method with the conventional Kramers- 
Kronig transform approach and find that the two methods agree over most of the 
frequency regime. For the same data set, the present approach is superior to Kramers- 
Kronig at high frequencies and can be made to perform at least comparable at low 
frequencies. 
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1. Introduction 

Passive micro-rlieology is a powerful experimental technique to probe the viscoelastic 
properties of hquids [H l2l [3] . The method can be used to perform in situ measurements 
of the hnear visco-elastic properties of complex hquids on samples that are many orders 
of magnitude smaller than those required in conventional rheology experiments. In a 
typical micro-rheology experiment, inert colloids are dispersed throughout the medium 
under study. The random (Brownian) displacements of the colloids are tracked optically, 
e.g. by video camera (~ kHz sampling rate) or using a quadrant photodiode (~ MHz 
sampling rate). In what follows, we focus on "real space" micro-rheology experiments, 
although it should be stressed that the key ideas behind the approach were first 
put forward by Mason and Weitz [1] in the context of diffusing-wave spectroscopy 
experiments. 

The starting point for all micro-rheology studies is the generalized Stokes-Einstein 
relation 



where D{s) is the Laplace transform of the time-dependent diffusion coefficient of 
a spherical particle with radius R and ri{s) is the Laplace transform of the time- 
dependent viscosity, which is related to the Laplace-transformed modulus G{s) by 
G{s) = sri{s). Eqn. [T]was proposed by Mason and Weitz [4J in the context of diffusing- 
wave spectroscopy and implemented for real-space (particle-tracking) micro-rheology 
by Gittes and Schnurr et al. [SI [1]. The approach was placed on a firm theoretical 
footing by Levine and Lubensky [6]. Excellent reviews of microrheology can be found 



In what follows, we assume that the conditions for the validity of Eqn. [T]are satisfied. 
As micro-rheology is by now a well-established and widely used technique, the focus 
of the present paper is not on the basic equations of micro-rheology, but rather on 
the most convenient way to translate experimental data on particle displacements into 
frequency-dependent, complex moduli. Again, several approaches to convert time series 
of particle displacements into complex moduli have been discussed in the literature 
(see, e.g. [HI E] HH 12]). However, the existing approaches suffer from one or more 
of the following drawbacks: 1) truncation errors in the transformation of the data for 
time series to frequency dependent moduli, 2) the need to fix the size of the data 
set (and hence the duration of the measurement) in advance, 3) the use of analytical 
approximations to the experimental data to facilitate the calculation of the visco-elastic 
moduli. 

The approach that we propose here has the advantage that it is fast and simple and 
it does not suffer from truncation errors. Analytical approximations are not needed to 
obtain an accurate high-frequency response. Statistical noise may affect the performance 
at low frequencies. In that case, the best results are achieved by using a properly 
normalized analytical approximation to the raw data. Moreover, due to the use of an 
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on-line data-reduction procedure, there is in practice no limit on the size of the data 
stream (and hence on the duration of the measurement). As a result, the emerging 
results for the complex moduli can be viewed "on line" during data acquisition. 

Moreover, the data-reduction procedure is effectively "loss free" - that means that 
one and the same data set can be used to probe the visco-elastic moduli at high and low 
frequencies. The choice of the low-frequency cut-off can be made during the experiments: 
by simply running longer, the low-frequency cut-off is decreased. However, we stress 
that the real advantage of the present method is not at low frequencies, but at high 
frequencies. This is significant because one of the key advantages of micro-rheology is 
precisely that it can access the visco-elastic moduli at frequencies that are too high to 
be probed in conventional rheology experiments. 

The basic idea behind our approach is to compute (during data acquisition) the 
time-correlation function of particle displacements and use a systematic coarse-graining 
procedure [13] that reduces the memory requirements for computing a correlation 
function from the full set of points to 0{hiN) points. Importantly, the method still 
makes use of all available data points and hence no physically meaningful information 
is lost in the coarse-graining process. 




646 




Figure 1. Video camera capture of a l/im diameter polystyrene microsphere in 
aqueous solution in an infrared (1064 nm) laser trap, and an example of tracking by 
image correlation. By tracking the Brownian motion of the particle, one can calculate 
the viscoelastic properties of the medium. 



Below, we describe our approach and show how, in combination with a variable- 
interval Fourier-Laplace transform method described in ref. [H], it can be used to 
compute complex moduli that agree with the results obtained by existing techniques 
(at least, in the region where the latter do not suffer from truncation errors). We first 
test our approach on "synthetic" data and then we apply it to experimental micro- 
rheology data. 
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2. Computation of complex moduli 

2.1. Relationship between time series and complex moduli 

We start from Eqn. [1] where the Laplace transform of the diffusion coefficient is given 
by 

POO 

D{s) = / dt e-''D{t) (2) 
Jo 

and 

D{t)= f dt' {vMv.{t')) (3) 





In Eqn [3] we have expressed the time-dependent diffusion coefficient as the integral of 
the "velocity autocorrelation function" (VACF) (fx(O)fx(t)). Note that the time scales 
in micro-rheology are such that one never measures the true velocity of the colloidal 
particle, rather one measures the (diffusive) displacement Ax(r) of such a particle 
in the shortest time interval r. In what follows, we interpret the "velocity" Vx{t) as 
{x{t+T) —x{t))/ T. With this definition - and the replacement of the integration in Eqn. [3] 
by a summation, we obtain the expression for the time- dependent diffusion coefficient 
that we will use below. Using the generalised Stokes- Einstein expression (Eqn. [1]), we 
can now relate the Laplace transform of the velocity autocorrelation function to the 
elastic moduli of the medium. Using eqns. [2] and [3] above, we can write 

/•oo /•£ 

D{s)= dte-^'l dt' {v,{0)v,{t')) 

where C{f(t)) denotes the Laplace transform of function f{t). As we shall see 
below, a similar expression can still be used when we consider a sequence of position 
measurements at discrete time intervals. Using eqn. [T] we can then write 

s QttRG{s) ^ ' 

If the colloidal particle is confined in an optical trap, then the force constant k of this 
trap appears as a correction to the elastic modulus: 

QtxRG{s) QttRG{s) + K (6) 

This expression is (of course) also correct in the absence of an embedding medium (i.e. 
with G{s) = 0), in which case 

D{s) = ^ (7) 
The mean-squared displacement of the particle is then 

((Ax)'(t^oo)) = 21ims^ (8) 



or 



kT 

{{Axnt)) = 2— (9) 



Real-time monitoring of complex moduli from micro-rheology 



5 



as expected from equi-partition, where £ ^ is the inverse Laplace transform. In the 
hmit K — 7- and G{s) = srj we recover the Einstein relation 

{{Ax)\t)) = 2Dt . (10) 

In general, our expression for G{s) is 

kT K 

^^'^ = QnRD{s) ~ 6^ ^^^^ 

Our key point is that it is advantageous to determine G{s) directly from the velocity 
auto-correlation function (or, more precisely, the discrete-time equivalent of this object). 
This correlation function is then coarse-grained online and D{s) follows directly from a 
single Fourier Laplace transform. 

Note that the conventional method to analyse micro-rheology data starts by 
computing the power spectrum of the particle displacements {\x{u)\'^). This power 
spectrum is related to imaginary part of the complex susceptibility a^u), that is related 
to the complex modulus G*{u) through 

G*(uj) = — 7^—-^ (12) 

where * denotes a complex quantity. The imaginary part of aiuS) is related to (|x(ci;)p) 
through 

Jil..\ ^ /l„r.M2 



«'V) = ^(kMr) (13) 

In order to obtain G{uj), we need to know both a' and a", a' can be obtained from a" 
via a Kramers-Kronig transform 

a [w) = -P / — -dC 

) 

cos{ut)dt I a"(C) sin(Ct)rf(C) (14) 



7r JO 

Note that the Kramers-Kroning route requires several transforms instead of one. 
Moreover, as is obvious from the integration limits of eqn. [TH the KK transform is 
sensitive to the high-frequency cut-off of the data. This tends to manifest itself as 
apparent unphysical behavior of the computed moduli at the limits of the frequency 
interval studied. Moreover, we shall show that the method based on the velocity 
correlation function allows a very convenient "on-line" reduction of the data, such that 
with very moderate storage {0{logN) - where N is the number of data points) - we 
retain the full information about the complex moduli over a frequency interval that 
extends from 27r/T to 27r/At, where T is the total time of the measurement and At is 
the interval between successive data points n + 1 and n. 

Below, we briefly describe how we perform this on-line coarse-graining of the 
incoming data stream. The method used is similar to the approach used to compute 
time correlation functions over very long time intervals in computer simulations [13] and 
for on-board data reduction in Dynamic Light Scattering experiments [15]. 
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2.2. Coarse- graining procedure 

Raw data, in the form of particle displacements, are continuously fed into a circular 
buffer of length N (e.g. = 100) in a cyclic fashion, such that once it is filled, the 
first data point (say point 1) is overwritten by the most recent data (say point N+1). 
The primary data are of the form 6x{nAt) = x{{n + l)At) — x{nAt), where x is the 
(calibrated) displacement of the particle that is being tracked and At the time interval 
between successive sampling points. The index n runs from 1 to Umax, where rimax is 
the total number of points sampled during the measurement. 

To achieve coarse graining, we sum the M most recently entered data points 
whenever n (mod M) = and enter the result into the next level 'coarse-grained' 
array that is also organised as a circular buffer. In this procedure, the size of the arrays 
for each level of coarse graining is preset and, as the number of coarse-graining levels 
grows only logarithmically with the total number of data points, overfilling the arrays 
is, in practice, impossible. The only important requirement is that the computation 
of the correlation function (see below) is faster than the data acquisition rate. In the 
example that we study, this is always the case. But even for data-acquisition speeds in 
the MHz regime, on-line processing of the data should not be a problem (for instance, 
by using a Graphics Processing Unit). If, for instance, we consider the case with A^=100 
points per array, and a coarse graining factor M=10, then the total storage required to 
accumulate the correlation function for a run of a billion points would be 2 x 8 x 100, 
including the "accumulator arrays" (see below). 



6x(i)(1) 5x(i)(2) 5x(i)(3) •■■ 5x(i)(M) 5x(i)(M+1) •■■ 6x(i)(N) 




5x(2)(1) 5x(2)(2) 5x(2)(3) •■■ 5x(2)(M) 5x(2)(M+1) •■■ 5x(2)(N) 




J 



5x(3)(1) 5x(3)(2) 5x(3)(3) ■■■ 5x(3)(M) 5x(3)(M+1) ■■■ 5x(3)(N) 



Figure 2. An example of data flow in a coarse graining scheme |T3] using an A^-point 
arrays and a coarse-graining factor of N/M. New data is cyclically introduced into the 
top (Level 1) array Sx'-^^t,), and averaged into Level 2 'coarse-grained' arrays with 
longer times Sx^^\t'j). These are subsequently averaged into an even higher Level 3 
fe(3)(i,) etc. 
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Data in the successive coarse-graining buffers are related such that 

M 

Sx^'^^\t) = Y.6x^'\nAt) , (15) 

n=l 

i.e. M data points from level i are summed and entered into the appropriate bin of the 
i + 1*^ buffer. 

In addition to the circular data buffers, we use an equal number of "accumulator" 
arrays where we store the information about the velocity auto-correlation function. 
Whenever a new data point is entered into an array of size N at coarse-graining level 
m, it is multiplied with itself and with all — 1 preceding points in the same buffer. 
The resulting N products are accumulated in a linear array of length N, such that the 
product of points N and N — n are added to bin n + 1. For every level of coarse graining 
there is a separate accumulator array. Computing the products and adding the 
resulting numbers to the appropriate elements of the accumulator array is the most time- 
consuming part of our data acquisition algorithm. However, even with MHz acquisition 
rates, modern processors (in particular CPU's) should easily be able to keep up with 
the data acquisition. 

The key assumption in the procedure described here is that the correlation functions 
that are measured show no high-frequency modulations at long times. This would clearly 
be incorrect if the particles would, for instance, perform a high-frequency, undamped 
periodic motion. Yet, for Brownian motion in a passive medium, this assumption is 
always satisfied. 

At the end of an experiment, we have an estimate of the desired "velocity" auto- 
correlation function. However, it is stored in different arrays - one for each level of 
coarse graining. Moreover, these arrays overlap. This is important because this allows 
us to check if the coarse graining introduces systematic errors at short times; for the 
first point, it always does. However, as we will show below, the effect of coarse-graining 
is minimal for subsequent points. 

The final step in our analysis is to construct a single correlation function out of the 
various coarse-grained arrays and apply a Fourier-Laplace transform to compute D{s). 
Each subsequent level of the correlation function has a time spacing that is M times 
larger than that at the previous level and extends for over a correspondingly longer 
time. By combining different levels of coarse graining, we are thus extending the low 
frequency limit of our micro-rheology data, while the high frequency limit is still set by 
the Nyquist condition. 

2.3. Direct conversion of coarse-grained time series into complex moduli 

In order to obtain G*{uj),we first need to compute the Fourier-Laplace transform of the 
VACF. Due to our coarse-graining procedure, we have values for the VACF at non- 
equidistant points in time. A convenient procedure to obtain the Laplace transform of 
such a data set was proposed by Evans et al. [E] . The method uses a linear interpolation 
to connect successive data point. Between data points, the slope of these linear segments 
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Figure 3. An example of a series of coarse-grained correlation functions, in this case 
VACFs calculated from simulated fluctuations of a Ifim sphere in an optical trap with 
trapping constant k ~ 10^^ pN/nm, diffusing with diffusion constant D = 10^^^ 
m^s~^. Correlations found from different levels are found to overlap. 



is (of course) constant, but at every data point, the slope may change discontinuously. 
The second derivative of this interpolated function is zero everywhere, except at the data 
points where the slope of the linear interpolation changes. There, the second derivative 
is a delta-function with an amplitude equal to the change in slope of the linear segments. 
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Figure 4. Transformation of an experimentally obtained correlation is found by 
considering a linear piecewise function which passes through all points and considering 
the functional form and magnitude of the first and second derivative. If C{t) is 
the original correlation function, C{t) and C{t) are the first and second derivatives 
respectively. 



C{ti 



'c{t,+,)-c{t,) c{t,)-c{t,^iy 



\ ti+l ~ ti 

ai5{t - ti) 



6{t- ti 



(16) 
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where Oj is a shorthand notation of the amphtude of the delta function ai t — ti. The 
Laplace transform of a series of delta functions is 

i 

for all ti^O. 

Next, consider the contribution of the point at t = 0. It is important to note that 
this first point of the VACF has an amplitude that is much larger than any of the points 
at t 7^ 0. For this reason, it is convenient to treat this point - and its Laplace transform 
- separately from that of the rest of the VACF. To this end, we spht the VACF into 
a delta-function part at t — and a remainder that has been defined such that it is 
continuous, and has a continuous first derivative, at t — 0. This is achieved as follows: 
we define Cs{t), the delta- function part at i = as 

Cs{t = 0) = {C{t = 0) - C{t = At)) At5{t) (17) 

where C{t — At)) is the value of the VACF after one time step. The remainder of the 
correlation function is the same as before, except that the value at i = is equal to 
the value at At. As a consequence, this function is continuous and has a continuous 

first derivative at t = 0. Hence, the second derivative of the linear interpolation of 
this function vanishes at t = 0. We first compute the Laplace transform of the part 
delta- function part. It is given by 

{C{t^O)-C{t^At)) ^^ 
2 

The factor 1/2 follows from the fact that the delta function is symmetric around t = 
when obtained as a limit of a delta function - we only integrate the part with t > 0. 
Note that integration of this part of the VACF yields a constant, real contribution to 
the diffusion coefficient, Dc- For confined particles D{t — > oo) = 0. However, due to 
statistical noise, this condition may not be satisfied. In that case, it can be enforced by 
adding a small correction to the first point of the VACF to ensure that D{t — > oo) = 0: 

C(t = 0) ^ C(t = 0) - 2 X ^ . 

Next, we consider Cvv{s), the Laplace transform of the remainder of the VACF. We 
use that fact that 

m = s'fis) - sF{0) - F(0) 

Hence, 

CMs) = ^+'^^^ + (18) 

s 

where the final zero on the right hand side follows because the first derivative of Cwit) 
vanishes at t — 0. It then follows that 

CMs) = ^^ + ^^^^ (19) 

s 
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This part of the Laplace transform of the VACF must be added to the constant 
contribution that resuhs from the transform of the delta function at t = 0. Combining 
fT9] and \17\ we get a complete Laplace transform of the VACF C(s), which is related 
to the diffusion constant D{s) using H] and thus the complex moduli via [TTl In what 
follows, we shall consider the Fourier-Laplace transform, i.e. we replace s by —tu. 

3. Test of method 
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Figure 5. The optical tweezer setup involved the use of a 1064 nm infrared laser, an 
AOD (acouto-optic deflector) unit and three stepping motors for moving the sample 
in the xy and z planes respectively. Visualisation was obtained using a CMOS camera. 



We tested the approach described above both on real and on "synthetic" micro- 
rheology data. The synthetic data were obtained numerically, using a kinetic Monte 
Carlo simulation of a particle in a harmonic trap. The objective of this test was to verify 
that the diffusion coefficient and the trap stiffness that were used as input parameters 
in the simulation, would be recovered in the subsequent data analysis of the VACF. 

Experimental data were obtained using an infrared optical tweezer setup that is 
described in ref. [16]. The laser beam was focused by a 63x water immersion objective 
with numerical aperture 0.90. The sample chamber consisted of a standard microscope 
slide, 100 fim SecureSeal imaging spacers (Grace Biolabs) and a coverslip. The chamber 
was inverted to ensure that the laser light was introduced from the cover-slip side. 
The sample was moved in the xy plane using a motorised stage, and an acousto-optic 
deflector (AOD) unit was used to accurately position the trap. Movement of the trap 
in the z direction was achieved using a separate motor moving the lens vertically. 
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The motion of microspheres was captured using a CMOS camera (Marhn F-131B, 
Alhed Vision Technologies) with a maximum frame rate of around 500 frames per second. 
This maximum rate was only achieved with the field of view restricted to one bead. To 
track particles, we used spatial correlation with an optimised kernel followed by a 2 
dimensional fit. Pixels were converted into distances using a graticule. 

The samples investigated were dilute (around 1 particle per 10 /im^) suspensions 
of 3/im diameter monodisperse silica micro-spheres (Microspheres GmbH) suspended in 
solutions of polyethylene glycol (PEG 8000) in MilliQ water at PEG concentrations of 
0, 20 and 40 % by weight. 

For the sake of comparison, we used bulk rheology measurements (Physica MCR 
benchtop rheometer (Anton Paar GmbH) ) to obtain independent estimates of G' and 
G" for these solutions. The viscosity data thus obtained were in good agreement with 
the literature data reported in ref . \171 . 

We used MATLAB for the data analysis and the Monte Carlo simulations. 

4. Results and discussion 




10^ 10^ 10^ 10* 10^ 10^ 

CO (Hz) 



Figure 6. Synthetic data for a Brownian particle moving in a harmonic trap. The 
elastic response of the trap results in an effective elastic modulus G'. In the same 
figure, we also show the results of the conventional analysis, using the Kramers-Kronig 
method. We note that the KK method suffers from systematic truncation errors at 
high frequencies. In contrast, the present method suffers from statistical noise at low 
frequencies. These problems can be addressed by replacing the (noisy) long-time part 
of the VACF by an analytical fit. See Figure [7] 



We start with the analysis of the "synthetic" data. We generated an artificial 



Real-time monitoring of complex moduli from micro-rheology 



12 



10 



10 



10" 



b 
b 



10" 



10" 



10" 



10" 





fit VACF 






/-VII 

fit VACF 




□ 




pts 


o 




pts 



10" 



10^ 10^ 
(o(Hz) 



10" 



Figure 7. Analysis of "synthetic" data for the same set of parameters as in Fig. [6] 
but with a data set that was 100 times larger (10^ points). The better statistics result 
in a significant reduction of the noise in the estimated visco-elastic moduli at both 
high and low frequencies. Note that, as in Fig.[6l G' is the "apparent" elastic modulus 
that is due to the stiffness of the harmonic trap. The drawn and dotted curves were 
obtained by fitting the VACF to a delta-function plus a single exponential. Clearly, 
such a fitting procedure further suppresses the effect of the statistical noise. 



trajectory of a Brownian particle in a harmonic trap, using the following parameters: 
diffusion constant D = 10^^^ m^s~^, temperature T = 298K and sampling rate 100 kHz. 
An artificial trap wsa imposed with k, = 10~^ pN/nm stiffness. As the harmonic trap 
was isotropic, moduli obtained from the X and Y displacements could be averaged. In 
most of our studies, the length of the data sets was 10"^ points (100s). To relate the 
diffusion constants to visco-elastic moduli, we chose a value of l/xm for the radius R of the 
microsphere. Using the Stokes-Einstein relation, a diffusion coefficient D = 10~^^ m^s~^ 
then corresponds to a viscosity r/ = 2.18 x 10""* Pas. The synthetic data correspond to 
a particle in a purely viscous liquid. In the case of a particle in an harmonic trap, the 
apparent elastic modulus G' follows from Eqn. [11] 

The moduli found are shown in Figure El In this figure, we have not corrected 
for the trap stiffness to obtain the true G' because we compare the effective G' with 
the value that follows from the known spring constant of the trap. The (apparent) real 
modulus is indeed found to agree with the predicted value Gjj-ap = 0.0531 Pa. G" follows 
the expected rjcu curve. 

To illustrate the difference between the present approach and the KK method, we 
have plotted the moduli as obtained by the Kramers- Kronig transformation method [5] 
in Figure [6l As can be seen from the figure, the KK approach becomes unreliable at high 
frequencies. This problem has also been reported by Schnurr et al. [1]. In contrast, the 
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present method appears robust in the frequency region where the KK method fails. This 
difference is important as the high-frequency samphng of the visco-elastic properties is 
one of the main strengths of micro-rheology: the use of a method that is robust at the 
highest frequencies is therefore important. 




CO (Hz) 



Figure 8. Calculation of G" from experimental micro-rheology data on pure water 
and on a 20% wt PEG8000 solution (see text). In this figure G' is not shown because, 
after correction for the trap stiffness, it does not differ significantly from zero. The 
straight lines are the predictions for G" based on the known low-frequency viscosities 
(see Ref. [IT]). 

At low frequencies, where G" ^ G', the visco-elastic moduli become noisy. The 
reason is that, whereas in the KK approach, a positive sign of G" is guaranteed even 
when G" ^ G', this constraint is not imposed for the present method. This problem is 
due to statistical noise and can easily be fixed by fitting the t > part of the VACF to a 
sum of exponentials with the constraint that the VACF (i.e. including the contribution 
at t = 0) integrates to zero (at least, for particles in a trap). Figure [7] shows how the 
use of such a constrained fit greatly improves the quality of the low-frequency moduli. 
Data fitting can be avoided altogether by simply collecting a larger data set. Here, we 
benefit from the fact that the storage required by our method scales only logarithmically 
with the size of the data set. As an illustration. Fig. [7] shows the results for G' and G" 
obtained from a set of 1 billion (synthetic) data points; analysing such a data set with 
the conventional (KK) algorithm would require at least 8 gigabytes of memory. 

Data sets were also taken for the suspensions of silica beads in 3 different 
concentrations of PEG 8000. In all three cases, we used the same strength of laser 
trapping as before. The bead diameters were 3/i for the and 20% measurements, l/xm 
for the 40% - the moduli should be independent of bead size. The curves shown in 
Figs. IS and |9] are the result of averaging moduli obtained using x and y displacement 
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Figure 9. Comparison of tlie complex moduli G' and G" of a 40% wt aqueous 
PEG8000 solution as obtained by classical rheology (closed symbols) and micro- 
rheology (open symbols). As this solution exhibits visco-eleastic behaviour, G' shows 
the (expected) quadratic increase with frequency, whereas G" agrees well with the zero- 
frequency viscosity at low frequencies (126.5 mPa s [17]) and levels off at near the point 
where G' peaks. As can be seen, the micro-rheology results are largely complementary 
to the conventional viscosity data. In regions of overlap, the two methods agree. Note 
that the micro-rheology data become very noisy for moduli less that 5 10~^ Pa. Longer 
runs (and/or fitting the VACF) would alleviate this problem. 



data. The error bars arise from polydispersity of the beads used. 

Interestingly, the viscosity of the most concentrated 40% wt aqueous PEG solution 
appears to decrease at high frequencies. This suggests that the stress autocorrelation 
function (that is related by a Green-Kubo relation to the viscosity) does not decay 
completely on the shortest time scales sampled in the experiments. This observation 
may be related to the non-Newtonian behaviour that has been observed in high 
molecular weight polyethylene oxide [18]. In the high frequency range, we also observe 
a corresponding increase in the elastic modulus with an intermediate ~ u"^ region, 
characteristic of a Maxwell fluid [T9] . 

At lower frequencies, the micro-rheology experiments yield data that are in good 
agreement with those that we obtained in bulk rheology measurements. Of course, the 
high-frequency range cannot be probed with bulk rheology. 

5. Conclusions 

In conclusion, we have developed a coarse-graining scheme that provides an efficient 
means to obtain complex moduli from micro-rheology measurements in real time. Even 
for very long measurements, the required memory usage of the method is very small. 
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This is important, as illustrated by recent work on the high bandwidth measurement 
of the VACF of a Brownian particle [20]. In this paper, it is stressed that the direct 
accumulation of the VACF is limited by the memory capacity of the data acquisition 
hardware. The approach produced here should present an easily implemented solution 
to this, and many other experimental studies related to micro-rheology. 
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